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Abstract 

Low energy cosmic ray antideuterons provide a unique low background channel for indirect detec- 
tion of dark matter. We compute the cosmic ray flux of antideuterons from hadronic annihilations 
of dark matter for various Standard Model final states and determine the mass reach of two fu- 
ture experiments (AMS-02 and GAPS) designed to greatly increase the sensitivity of antideuteron 
detection over current bounds. We consider generic models of scalar, fermion, and massive vector 
bosons as thermal dark matter, describe their basic features relevant to direct and indirect detec- 
tion, and discuss the implications of direct detection bounds on models of dark matter as a thermal 
relic. We also consider specific dark matter candidates and assess their potential for detection via 
antideuterons from their hadronic annihilation channels. Since the dark matter mass reach of the 
GAPS experiment can be well above 100 GeV, we find that antideuterons can be a good indirect 
detection channel for a variety of thermal relic electroweak scale dark matter candidates, even when 
the rate for direct detection is highly suppressed. 
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I. INTRODUCTION TO ANTIDEUTERONS 



Antimatter is rare in our Universe, making it potentially a good low-background channel 
in which to search for evidence that dark matter annihilates to Standard Model particles. 
Due to color multiplicity, it is reasonable to expect a sizable fraction of annihilations into 
hadronic final states-either directly or through subsequent decay. Dark matter annihilation 
to hadronic final states can account for the antiproton {p) cosmic ray spectrum, which 
also provides constraints. This component has been calculated for various models of dark 



- |l3| constrains the dark matter's mass 



matter And antiproton cosmic ray data 

and annihilation rate for models with antiprotons in the final state. For instance, for dark 
matter with a mass tjidm < 50 GeV, {(y\v\)therm ~ 1 pb (cross section required by thermal 
relic density), and annihilation to bb, one finds that the dark matter contribution to the p 
spectrum can exceed experimental bounds jf], [s, [l^. On the other hand, for dark matter 
masses of moM > 50 GeV and an annihilation cross section {(j\v\)ann ~ 1 pb, the primary 
contribution to the antiproton cosmic ray flux of dark matter annihilation will generally be 
smaller than the secondary contribution from astrophysical processes. Therefore, only if the 
antiproton cosmic ray flux is known to a high degree of precision would it be possible to 
disentangle a dark matter signal from the astrophysical background. 

Compared to antiprotons, antideuteron cosmic rays are poised to be a more sensitive 



indirect probe of dark matter with a hadronic annihilation channel 15|-19|. The threshold 
energy required to create an antideuteron from cosmic ray proton collisions with interstellar 
H and He is Efh = 17mp ^ (for antiprotons it is Eth = 7mp). Furthermore, the low value 
of the antideuteron binding energy {B = 2.2MeV), removes its ability to slow down as it 
propagates in the galactic disk, leading to a low background for low energy (T5/„ < IGeV) 



antideuterons having an astrophysical origin |15l . l20l . l2lj . The fact that antideuterons can 



act as a more sensitive probe of the hadronic final states of dark matter annihilation was 
first investigated in 15^ where an antideuteron signal originating from supersymmetric dark 
matter was investigated. 

Future experiments, two energy bands of the Alpha Magnetic Spectrometer Experiment 



^ A pp collision which makes an antideuteron must have s > SGm^. In the rest frame of a Hydrogen nucleus 
in the galactic disk s = {irip + Ep)'^ — fc^, where Ep and kp are the energy and momentum of an impinging 
cosmic ray proton. This implies that Ep > I7mp in order for this process to occur. 
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AMS-02) 22|, |23| and the three ^ phases of the General Antiparticle Spectrometer (GAPS) 



25| , should greatly increase the sensitivity of searches for cosmic ray antideuterons over 



;he current bound on the antideuteron component of the cosmic ray Spectrum set by BESS 



12| . The current bound is: 

$5 < 0.95 X 10-^[m2 s sr GeV]-^ 0.17 < T/n < 1.15 (GeV/n) BESS (1) 

The proposed sensitivites of the future AMS-02 and GAPS experiments are: 

$5 = 2.25 X 10"^ [m^ s sr GeV]-^ 0.2 < T/n < 0.8 (GeV/n) AMS-02, (2) 

$5 = 2.25 X 10"^[m2 s sr GeV]"^ 2.2 < T/n < 4.2 (GeV/n) AMS-02, (3) 

$5 = 1.5 X 10"^[m2 s sr GeV]-^ 0.1 < T/n < 0.2 (GeV/n) GAPS(LDB) (4) 

$5 = 3.0 X 10-^[m'^ s sr GeV]"^ 0.05 < T/n < 0.25 (GeV/n) GAPS(ULDB) (5) 

$5 ~ 2.6 X 10-^[m2 s sr GeV]-^ 0.1 < T/n < 0.4 (GeV/n) GAPS(SAT). (6) 

For higher values of antideuteron cosmic ray flux than those quoted above, antideuteron 
detection is expected. Since the contribution to $^ from astrophysical background is 
~ 5 X 10^^ [m^ s sr GeV]^-^, these experiments (with the exception of a satellite 
mission) will be essentially background free probes of an anomalous component of the an- 
tideuteron spectrum. Any detected signal must come from some source that contributes 
significantly to the low energy antideuteron spectrum. Dark matter can generate such a 
signal. 

If dark matter annihilates, this process essentially occurs in the Earth's rest frame because 
the halo and Earth frame of reference are approximately the same on relativistic scales. 
Therefore dark matter annihilation can act as a source of antideuterons that populates the 
low energy antideuteron spectrum. If the dark matter annihilates at a large enough rate, 
then this signal can be detected at the AMS-02 and GAPS experiments. Furthermore, as 
was investigated in |l6|, if the Dark matter annihilates directly to colored states (e.g. qq), 
then the final state spectrum of antideuterons tends to be peaked at low energies due to 
the fact that hadronization for the final state occurs in the rest frame of the halo, but when 
the hadronization takes place in a boosted reference frame, for instance if the final states 



^ The three proposals are the Long Duration Balloon(LDB), the Ultra-Long Duration Balloon(ULDB), 
and a satellite (SAT) mission. 
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are boosted color singlets (e.g. W~^W~), the antideuteron spectrum tends to be peaked 
away from low energy. If the color singlet final states are produced at threshold however, 
the antideuteron spectrum will be peaked at low energies. So, not only do antideuterons 
provide an essentially background free channel in which to indirectly probe dark matter 
annihilations, but the detection of a signal is preferred by some annihilation modes compared 
to others. 

It is common in many Beyond the Standard Model (BSM) physics scenarios for the 
dark matter to have a significant annihilation channel with hadrons in the final state. An- 



tideuteron signals from supersymmetric neutralinos were investigated in 



15 



3; Q 



in- 



cluded an analysis of the antideuteron signal from Universal and Warped Extra dimensional 



dark matter candidates. In [18| the antideuteron signal was given as a function of dark mat- 
ter decay rate. In heavy doublet dark matter with an enhanced present day annihilation 
cross section was investigated. Independent of a particular BSM physics model, the exper- 
imental reach of the GAPS experiment as a function of the present day annihilation cross 
section and mass for the annihilation channels of WW and bb final states was investigated 
in Q- 

The first goal of this paper is to present a model independent assessment of the dark mat- 
ter mass reach of the AMS-02 and GAPS experiments for different annihilation channels. 
In addition to the bb and W^W" final states, which have appeared in previous analysis, 
we also assess the reach for tt, hh, and gg final states as a function of dark matter mass. 
Furthermore, for our injection spectrum computation, we utilize the non-fact orized "coales- 
cence scheme" recently advocated by 261], which more accurately characterizes antideuteron 
production for larger kinetic energies and has significant effects even for the low energies 
of the GAPS and AMS-02 experiments. In Section [TTl we describe our calculation of the 
antideuteron cosmic ray flux. In Section [Till we present our results as a number of expected 
events as a function of dark matter mass for a given final state annihilation mode. In Section 
IIVI we describe the apparent conflict between a thermal relic annihilation cross section and 
bounds on the spin-independent elastic scattering cross-sections derived from direct detec- 
tion experiments. We then describe the basic mechanisms utilized by many models that 
evade this conflict. 

The second goal of this paper is to understand the status of some well-motivated spe- 
cific dark matter candidates considering both direct detection and indirect detection via 
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antideuteron cosmic rays. In Section |V] we apply our model independent results to a small 
sample of dark matter candidates, describe each model's essential features by analyzing the 
leading operators which govern present day annihilation and elastic scattering off of nuclei, 
and determine the extent to which they can be probed via antideuteron searches. Finally, 
in Section IVII we conclude. Appendix A contains reference material for understanding the 
kinematic behavior of annihilation cross-sections and elastic scattering cross-sections for 
general DM candidates based on operator analysis. Appendix B contains a brief description 
of an inelastic dark matter candidate analyzed in SectioifVl which is compatible with usually 
stringent solar neutrino bound while explaining DAMA signal. 



II. ANTIDEUTERON COSMIC RAY FLUX 



The antideuteron cosmic ray flux at Earth can be decomposed into primary, secondary, 
and tertiary components. Standard astrophysical sources produce only secondary and ter- 
tiary components; these originate from astrophysical scattering o 



protons off the interstellar 

1 n 

Hydrogen and Helium and have been thoroughly studied in [20|, |2l[ and updated in [17]. A 
primary component can potentially come from the hadronic annihilation channels of dark 
matter. Here we focus on contributions from various types of dark matter annihilation. 

Only three properties of dark matter determine its contribution to the antideuteron cos- 
mic ray spectrum, its mass ( rriijM ), annihilation cross section ( {cr\v\)ann ), and fraction 
of annihilation into a particular final state. If there are multiple annihilation modes, their 
contribution just adds linearly to a good approximation. Connecting these dark matter 
properties to the antideuteron flux at earth requires three steps. First, one must calculate 
the injection spectrum of antideuterons from dark matter annihilation. Second, one must 
propagate the antideuterons according to a consistent model of cosmic ray propagation. 
Third, one must account for the effects of solar modulation. In the following subsections we 
review the calculations involved in these three steps. 



A. Injection Spectrum 

In order to calculate the injection spectrum of antideuterons from a particular dark 
matter model, one needs to know the dark matter mass and annihilation cross section into 
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various final states. The dark matter mass sets the energy scale of the final state, and 
the composition of the final state determines how that the system hadronizes into (anti-) 
baryons. This process can be modeled by numerical Monte Carlo. We choose PYTHIA 
6.400 j^^] for our computation which were run on the Odyssey cluster supported by the 
FAS Research Group. Then one must make a prescription as to how a p and a n form 
an antideuteron. Antideuteron production, whether from p + H^D + X or dark matter 
annihilation requires such a prescription. One prescription is the so-called "Coalescence 
Model" 20| . According to this model in the n rest-frame of a n — p system, if the p has 



kinetic energy less than the antideuteron binding energy, B, then one assumes the two bind 
to form an antideuteron. This crude model provides a good order of magnitude estimate 
for the formation of antideuterons. Said more precisely, if KEp B {\kp\ ~ {2mpB)^) then 
antideuteron formation is assumed. So the condition that a p and in a hadronic final 
state form an antideuteron is: \kn — kp\ < {2mpB)^ ~ ~ 70 MeV. po is referred to as 
the "coalescence momentum." We let this provide us with the conceptual understanding, 



but fix Po to reproduce t 
hadronic Z-boson decay 



le^ALEPH measurement of (5.9 ± 1.9) x 10 ^ antideuterons per 
. Typically, when using the "coalescence model" one assumes 



that the distribution of n and p with energy Ep^n is spherically symmetric and factorizes. 
sospin symmetry is also usually assumed in the analysis. We will follow the analysis of 



26| and compute the coalescence condition on an event by event basis. This method avoids 



the assumption of spherical symmetry and does not neglect phase space correlations in the 
production of a. n/ p pair. Just as in [26|, we find that a value of po = IGOMeV is required 
to match the ALEPH measurement; we use this value in all calculations that follow. 

We have computed the injection spectra for antideuterons for specific final states. In 
Fig. [1] we have plotted the injection spectrum for antideuterons for tt, bb, h^h^ , gg, and 
W^W~ final state annihilation channels for various masses. One can see a clear distinction 
between colored and color neutral final states for larger values of the dark matter mass. 
The colored final states yield a spectrum that remains peaked at low energies as the mass 
increases. However the color neutral final states yield spectra that are peaked at higher 
energies as the mass increases. This is due to the fact that hadronization takes place in a 
boosted reference frame for energetic color neutral final states. But notice that when the 
color neutral states are produced near threshold, the fact that they are nearly at rest implies 
hadronization takes place in the lab frame resulting in a spectrum not too different than the 
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colored final states. Note also that the tt spectrum is generally larger than the bb spectrum 
for the same dark matter mass because t ^ b + before hadronizing. Then, in addition 
to the bb component of the tt final state, when the decays, and 2/3 of its decay products 
are hadrons, the total final state will contain more hadrons, and thus more antideuterons, 
in general. A related feature of the top injection spectrum is that it generally contains a 
sightly higher energy antideuteron component due to the fact that as the dark matter mass 
increases, the Ws carry more of the initial top's kinetic energy and tend to populate the 
higher energy component of the antideuteron spectrum. A similar disparity exist between 
then W^W~ and h^h^ spectra but here it is simply that for a light higgs — t- bb 100% 
while — )■ hadrons only 2/3 of the time. So, for a given dark matter mass the light higgs 
final states have a higher antideuteron yield than do the Ws. There is also a small effect 
from the fact that the higgs are heavier than the Ws and thus are never quite as boosted 
as the Ws. The distinction between the two spectra becomes more dramatic when the higgs 
mass is significantly larger than the W-mass. For heavier higgs {mhiggs > 140 GeV), it 
decays dominantly to WW or ZZ and one finds that the spectrum resembles the W final 
state spectrum as the mass increases. In order to simplify our analysis, in Fig. [1] as well 
as for the remainder of the paper, we fix mhiggs = 115 GeV. We have not given explicit 
injection spectra for ZZ and light quark final states because these are very similar in shape 
and magnitude to the W^W~ and bb injection spectra, respectively. 

Finally, we comment that the shape of these injection spectra are fiatter compared to 
the injection spectra one derives from using the factorized "coalescence model" due to the 
fact that higher energy antideuteron production is not suppressed using an event by event 



calculation 



26|. 



B. Dark Matter Profile 



The shape of the dark matter halo profile can be parameterized as an NFW 291, Einasto 

n n 

[30|, or Isothermal [31] profile. We will restrict our attention to the Einasto profile due to 
the fact that it provides the most conservative signal while still being preferred by numerical 
simulations. 



pEin[r) = PqBXP 



-2 
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FIG. 1: The antideuteron injection spectrum as a function of Kinetic Energy, T, for 
DM annihilation to W~^W~ , hh {mhiggs = 115 GeV), it, bb, and gg final states. For 
each final state we plot the injection spectrum for dark matter masses of: rriuM = 
100 GeV (blue/solid), 200 GeV (green/dashed), 300 GeV(red/dotted), 400 GeV (black/dashed - dotted), 
500 GeV(black/solid),600 GeV(blue/solid), 700 GeV(green/dashed), 800 GeV (red/dotted) as com- 
puted using PYTHIA 6.400 and a coalescence momentum pQ = 160 MeV. However, for hh and 
ti final states the lowest dark matter mass plotted is = 150 GeV(blue/solid) and 

mriDM = 200 GeV(blue/solid), respectively. 
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Model 


6 


Ko (kpcVMyr) 


L(kpc) 


Vc (km/s) 


MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.70 


0.0112 


4 


12 


MAX 


0.46 


0.0765 


15 
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TABLE I: Propagation models 



with Pq = 0.4 GeV/cm^, = 20 kpc, and r© = 8.5 kpc. 



C. Propagation 

In order to connect the injection of antideuterons in the galactic halo to the flux of 
antideuterons near the Sun, one must propagate the antideuterons in the galactic halo. For 
propagation, our computations use a two dimensional diffusion model and our discussion in 



this section follows the notation of 



32 



34| (for reviews see [35|, |36| ) . In order to parameterize 
the uncertainty in propagation, we compute CR propagation with MIN, MED, and MAX 
propagation parameters as listed in Table [land given in jsj. 

The diffusion equation for charged cosmic rays, neglecting energy losses, is given by 

j^iP{r,z,E) = Q{r,z,E)-2hSiz)TanniE)inH + ^^nHe)Hr,z,E) 

where ilj{r, z, E) is the number density of a cosmic ray species. Steady state solutions 
are given by setting ^'^^ir, z^E) = 0. Also, the region over which this equation is solved 
is assumed to be a cylinder with height L, which can vary as in Table [H and a radius 
fixed at i? = 20kpc. A conservative estimate of the CR flux can be made by assuming 
that ip{R, L, E) = 0, at the boundary. A second boundary condition is given by assuming 
continuity through the z = plane. Solutions are expressed in terms of a Bessel Series 
Expansion. 

It is common to write fluxes in terms of particle kinetic energy (T) rather than energy 
(E), so we do this in what follows. The function Q{r, z, T) depends on the injection spectrum 
as 
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The Hydrogen and Helium number density in the galactic disk are: uh = 1/cm^ and nne = 
.07 uh, respectively. TanniT) describes the depletion of the antideuteron population due 
to dissociation after scattering from interstellar Hydrogen and Helium^, TanniT) = a{T)v 
where 

661(1 + 0.0115T-°-^^4 - 0.984T° °i^Mmb T < 15.5 GeV 
a{T)^{ ^ ' (10) 

36T-°-5mb, T > 15.5 GeV 

K{T) = K,{R{T))\ i?(T) = , E{T) = T + m (11) 

with T in units of GeV. Kq sets the overall scale for charged particle diffusion and Vq 
accounts for a "galactic wind" that pulls charged particles out of the disk. Kq and Vc for 
different sets of propagation parameters are hsted in Table I. The solution for the number 
density of CR's is given by: 

V'(it:o,o,r) = J]7v,(o,r)jo(^eO (12) 



where 



with 



ViiT) = 2 X sinh (|(L - z'))q,{z')dz' , (14) 



i^2(T)+^^J ' ^^ = 2/ir„„„ + yc + i^(T)coth(S,L/2) (16) 
where are the zeros of Jq. The cosmic ray flux at the Sun is given by 

^Rq,0,T) = ^^{RqAT). (17) 



D. Solar Modulation 

Finally, to relate the flux at the Sun's position to the flux at Earth, we must take into 
account effect of solar modulation. Solar modulation reduces the flux of low energy cosmic 



^ Note that F, as it is used here, does not stand for a decay width. 
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ray species. This final step is of particular importance since we are interested in the low- 
energy part of the antideuteron spectrum, where solar modulation effects are most important. 
We take into account solar modulation using the Gleeson Axford force field model |37| : 



We)= 2mT + T? ^^^^' T = T^ + e<P^. (18) 



2mT0 + T2 
2mT + J 

We will take (pp = 500MV in what follows. 
E. Backgrounds 

When considering whether or not a predicted signal from a primary source is large enough 
to be significant or not, it is important to understand the astrophysical sources which gener- 
ate the background. The background computation of secondary and tertiary antideuterons 
due to protons colliding with interstellar Hydrogen and Helium were originally carried out 



m 



20 



2l| . For the analysis that follows, we use the expected antideuteron flux at the top 



of the atmosphere for "MED" propagation parameters as computed in 17l |. 



III. EXPERIMENTAL REACH FOR CERTAIN FINAL STATES 

We now propagate the various injection spectra according to the procedure outlined 
above, assess the sensitivity of antideuteron search experiments, and comment on uncer- 
tainties in the results. 

A. Results 

Modulated antideuteron fluxes and experimental sensitivities are presented in Fig. [2J The 
fluxes are plotted for a fixed dark matter annihilation cross-section of {cr\v\)ann = 1 pb and 
various dark matter masses and annihilation modes and "MED" propagation parameters. 



The predicted background antideuteron flux as calculated in [17| is also plotted. We shall 
first discuss the features of the flux originating from dark matter annihilations (signal) and 
then consider the effects of the background when we estimate the experimental reach of the 
AMS-02 and GAPS experiments. 

Here we see that for each final state there exists some mass range for which AMS-02 and 
GAPS are sensitive to that particular dark matter annihilation mode. Let us emphasize that 
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FIG. 2: The antideuteron reach of the AMS-02 (black dotted) and GAPS (ULDB) (black upper- 
Dashed) and GAPS (SAT) (black lower-dashed) experiment for dark matter annihilation to 
W~^W~, hh {rrihiggs = 115 GeV), it, hh, and gg final states. In each case the dark matter present 
day annihilation is set to (cr|v|) = 1 pb. The fluxes are plotted for different dark matter masses 
from m_DM = 100 — 800 GeV as in Fig. [TJ The flux decreases as the mass increases. Also, the 
astrophysical background is given by the solid line. Propagation is done with "MED" parameters. 
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FIG. 3: The expected number of primary antideuterons for the AMS-02 (high) (blue/sohd), 
AMS-02 (low) (green/dashed, )GAPS (LDB) (red/dotted), GAPS (ULDB) (black/dot-dashed) 
and GAPS (SAT) (pink/soUd/upper hne) experiments for dark matter annihilation to W^W~, 
hh {rrihiggs = 115 GeV), it, bb, and gg final states. {cr\v\)ann = 1 pb- Propagation is done with 
"MED" parameters. 
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for even the largest antideuteron signal listed here, the antiproton spectrum is predicted to 
remain consistent with current BESS and PAMELA data, this is true even if one uses the 
"MAX" set of propagation parameters for antiprotons. Assuming that the sensitivity lines in 
Fig. |2] correspond to the fluxes required for a detection of one antideuteron, we convert each 
flux in Fig. [2] into a number of expected antideuteron events from primary sources {N^^"^) 
at each of the experiments (AMS-02 (high), AMS-02 (low), GAPS (LDB), GAPS (ULDB), 
GAPS (SAT)) and plot the results in Fig. [31 From Fig. |3l it is clear that the antideuteron 
signal for dark matter annihilation to colored final states does not fall off as rapidly as for 
color neutral states. This reflects the fact that hadronization and fragmentation takes place 
in a boosted frame of reference for higher energy color neutral final states. From this plot we 
can also see that the predicted number of antideuteron events can be 10 for some scenarios 
and even 100 for a GAPS level sensitivity on a satellite. However, for larger masses the 
number falls quite rapidly. 

Notice, again, that the shape of the flux differs from the flux one would derive using the 
factorized "coalescence model" due to the fact that the injection spectra predicted from the 
two models differs for larger antideuteron kinetic energies. In this case the flux is sustained 
at higher energies. This is important when considering experiments like AMS-02 which are 
sensitive to antideuterons with higher kinetic energies. 

We parameterize the background as recently given in jl^] which builds on the work of 



20 



2l| and have inferred from their results and the quoted experimental sensitivities, eq. (^^, 



that the expected number of background events at each of the experiments is: 

AMS - 02(high) : b = 0.3, AMS - 02(low) : b = 0.1, 

GAPS(LDB) : b = 0.04, GAPS(ULDB) : b = 0.2, GAPS(SAT) : 6 = 3.5 (19) 

In order to characterize the mass reach of antideuteron search experiments we use a 



statistical analysis similar to that of [16|], where we assume a Poisson distribution for the 
detection of antideuteron events. In Table [III we present, for various annihilation modes, 
the dark matter mass for which the number of expected antideuteron events equals Ncrit, 
where Ncru is determined as follows. Ncru is the lowest value of satisfying the following 
inequahties: 



Af-l N-1 



^P{n,b) > 0.997 (3a), or ^ P(n, 6) > 0.9999994 {5cr) (20) 



n=0 n=0 
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160 


2 


GAPS (SAT) (3a) 


500 


700 


500 


900 


240 


10 


AMS-02 high {5a) 


30 


< 'mt 


< 


60 


< my/ 


6 


AMS-02 low {5a) 


70 


< mt 


< rnh 


140 


< mw 


4 


GAPS (LDB) {5a) 


75 


< m 


< irih 


150 


< mw 


3 


GAPS (ULDB) {5a) 


150 


220 


150 


300 


120 


5 


GAPS (SAT) (5a) 


360 


550 


300 


670 


200 
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TABLE IL Mass reach for various experiments (all masses in GeV units). Masses quoted denote the 
dark matter mass for which the signal -|- background flux predicts a number of events equal to Ncrit- 
See text for a definition of Ncrit- For each annihilation mode, 100% branching, {a\v\)ann = 1 pb, 
MED propagation parameters, and coalescence momentum of po = 160 MeV is assumed. 

where P{N,b) = {h^ e~^) / N\ and h is the expected number of background events. In other 
words the probability of detecting the number of antideuterons listed in the Ncrit column of 
Table [III for each experiment is: 

Aa(iV > iV„,t) < .003, or P5a(iV > iV,,,i) < 0.0000006 (21) 

In these detection of Ncrit antideuteron events implies that one is seeing an exotic 

contribution to antideuteron cosmic rays at a confidence level > 95% and > 99.9999%, 
respectively. The masses in Table Hi] are, given an annihilation mode, the largest mass for 
which the antideuteron flux arising from the sum of the dark matter component and the 
background yields an expected number of events equal to Ncrit- The number of primary 
events expected from a dark matter source can be inferred from Fig. [31 This means that 
if dark matter is annihilating to ti final states with a mass m^iM < 400 GeV or hh with 
rriDM < 250 GeV then one should expect to see a significant source of primary antideuterons 
at GAPS (ULDB). Compared to the fact that the antiproton measurements bound dark 
matter masses of 50 GeV [16] with the same annihilation channel, GAPS(ULDB) will be 
able to probe dark matter masses which are greater by an order of magnitude. 
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B. Uncertainties 



Theoretical uncertainties arise at various points in the computation that produced the 
results in Fig. |2] and Fig. [31 They are uncertainty in the coalescence momentum (po) that 
determines the probability that an antideuteron forms during hadronization, the cosmic 
ray propagation parameters, the shape of the dark matter halo, and solar modulation. A 
discussion of these uncertainties can also be found in [16] and [l7|. 

The predicted antideuteron flux is most sensitive to the choice of "coalescence momen- 
tum" (po) used to calculate the injection spectrum of antideuterons from dark matter an- 
nihilation. Because the "coalescence model" keeps only events within a 3-d volume of 
phase space set by the parameter po, the overall flux of antideuterons depends on po as: 
~ Po- The choice and uncertainty of po is therefore very important. Recently, 26 1 
has argued an event by event computation of the coalescence condition is essential for a 
proper computation of higher energy injection spectra. After matching to ALEPH's mea- 
surement of (5.9 ± 1.9) X 10~^ antideuterons per hadronic Z-boson decay 28||, the resulting 
Po = 160 ± 17 MeV, translating into an uncertainty in the flux of 90%. This is the uncer- 
tainty estimate for simulating dark matter annihilation for a dark matter mass of 45 GeV. 



Based on the similar values for po deduced at lower energies [21|, |38| compared to those 
inferred from the ALEPH data, we expect that this uncertaintly in po should not increase 
significantly as the center of mass energy of the collision increases to 1 TeV. 

The next significant uncertainty in the antideuteron flux comes from the uncertainty in 
the parameters of cosmic ray propagation listed in Table [B Roughly, these amount to a 
variation in the flux by an order magnitude above and below the allowed flux predicted 
from using "MED" propagation parameters. 

The dark matter halo proflle is also uncertain with NEW (Navarro-Erenck- White) , 
Einasto, and Isothermal proflles as standard choices. When the dark matter density, p, 
becomes large, dark matter annihilation is enhanced as p^. Since the flux of antideuterons 
is sensitive to the halo shape at the center of the galaxy, the flux can depend strongly on 
the choice of halo proflle. In |l3], it was shown that there is a correlation between the 
uncertainty due to the halo proflle and the choice of propagation model. Eor MED and 
MIN propagation parameters the uncertainty in the low energy region of antideuteron flux 
is (9(5%) or less, while for MAX propagation parameters this uncertainty can be as large as 
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0(50%). 

Overall the uncertainty in the theoretical prediction for the comic ray antideuteron flux 
is dominated by the uncertainties of cosmic ray propagation. The actual flux could be a 
factor of 10 larger or smaller than the average values used in this paper. 



IV. GENERAL BOUNDS /FEATURES OF THERMAL DARK MATTER 

In the previous section, we have essentially done a 'model-blind' study of antideuterons 
as a probe of thermal WIMP DM because we have assumed (cr|f = 1 pb and varied the 
DM annihilation mode and mass. In the remainder of this work we will apply this general 
analysis to specific DM candidates. Much existing work on antideuteron DM searches focus 
on the reach for particular, well-motivated models like SUSY LSP ISMlTj]: we shall consider 
these and other models in the following sections. 

Before addressing specific models, in this section we will start with a very broad view and 



examine the conditions under which general WIMP DM can simultaneously yield t 



le right 



thermal relic density and satisfy current direct detection (DiDt) bounds. See [39144 1| for 
related discussions. The relevant DM properties are the DM quantum numbers (e.g. charge 
and spin) and the type of interaction it has with the Standard Model. DM properties can 
be organized by considering general DM bilinear operators and the Standard Model bilinear 
operators to which DM couples. We leave the detailed discussion and summary of these 
operator properties to Appendix A and will reference the results listed there as needed 
during the following sections. The important point that one gains from such an analysis 
is that a variety of thermal WIMP DM models can have very suppressed DiDt rate while 
still providing a good signal for indirect detection (IdDt) via antideuterons. The general 
discussion here and in Appendix A also provide a more systematic way of understanding 
the specific models we will consider in the next section. 



In general, there are four important cross-sections for a thermal dark matter candidate. 
The first is its annihilation cross section at early times {cf\v\)therm- This determines the 
thermal relic abundance of dark matter. For this work we are assuming that {(y\v\) therm = 
1 pb which leads to Q^Mh"^ = 0.12. The second is the present day annihilation cross- 
section {cr\v\)ann- This is important for indirect detection. The third and fourth are spin- 
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independent (SI) and spin- dependent (SD) scattering cross section off nucleons: asi and 
(JsD, respectively. These are important for direct detection. 

In this section, we first discuss the general correlation between {<j\v\) therm and a si and how 
current constraints on asi from Direct Detection (DiDt) experiments imply the inequality: 

^'^l^l^^^'^"" > 10^ (22) 

O'SI 

for a thermal relic of EW scale mass. First, we elaborate on why a thermal relic must 
satisfy this inequality. Then we discuss the essential features allowing models to satisfy 
this inequality. In particular we discuss how the physics determining {<j\v\) therm can be 
significantly different than asi- 



A. Thermal Relic Density versus DiDt Bounds 

Assuming that Q^Mh"^ = 0.12 arises because dark matter is a thermal relic implies that 



{(^\v\) therm = 1 pb at frcczcout. The CDMSII j42| and XENONIOO [43] DiDt experiments 
bound the spin- independent (SI) scattering cross section of DM on target nucleon of roughly 
asi ^ 10^'' pb for a typical WIMP in the mass range ~ 10 — 10'^ GeV. In general these 
two cross sections can be correlated via crossing symmetry of the Feynman diagram that 
controls dark matter annihilation. Consider, for example, a scalar dark matter candidate 
{x and its antiparticle) which annihilates to quarks, leptons, W/Z with 'unbiased' universal 
couplings, has DM mass m^, mediator mass M, width Tm, and mediator couplings to 
dark matter and Standard Model state (71,(72 (for scalar DM gi is of mass dimension 1), 
respectively. Then we get an effective Fermi coupling for the related operator x^XQQy 

r = ^ (o-i) 

[(4m2 -M2)2 + r2^M2]V2- ^ ^ 

Let the total annihilation cross section be {a\v\)therm = 1 pb- In order to relate this to DiDt, 
we focus on processes involving the up quark. With a color factor and universal couplings, 
it is reasonable to assume that the annihilation fraction to u quarks is 10%. 



A'elMP 



where after spin sum, |Mp ^ SG'^m'^. So, 
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Here we get a required relation 



{9192 



i2 



3 X (26) 



47r[(4m2 - M^f + Tl^M^] 
Crossing the Feynman diagram we get an associated process for DiDt (here we need to sum 
over all light quarks and the heavy quark contribution via a gluon loop). Take scattering 
off the proton as an example. The cross section is ^ 

_ 1 ml (g^g^f f rup ^ mp2_ \ 

^ ^ ' ' \q=u,d,s ^ q=c,b,t ^ / 

_ 1 {9i92? 

where f^Q, f^^ are proportional to gluon and quark matrix element in the nucleon. To get 
the second line we have used: ^ nip, with the measured values of: f^Q, f^^: /-^^ = 
0.020 ± 0.004, = 0.026 ± 0.005, /^^ = 0.118 ± 0.062, /^^ = 1 - J^^as fxq which is about 
0.84. 

Meanwhile the current bounds from XENONIOO and CDMS-II imply for ~ lOOGeV 

-4%^ ^ 3 X 10--cm^ (28) 
However, thermal relic density requirements (eq.f l26p ) imply 

asi = a^p = j^, ^3 X lO^^^cm^ (29) 

For EW scale m-^ and M, this implies cr^p ~ 10~^^cm^, which is about 300 times above 
current the DiDt bound. This bound weakens if > ITeV with a heavier mediator. 
The above argument suggests that one would naively expect 

lOl (30) 



{o-\v\)therm . „5 



So at least one of the assumptions we have made above must not hold in order to sat- 
isfy eq. f l2^ . for any realistic dark matter candidate. There are essentially two ways that 
eq. can be satisfied, by either having {a\v\) therm = 1 pb as a result of some kinematic 
enhancement or having some suppression for asi- 



^ The quantity DiDt experiments bound is the DM-nucleus scattering cross section per nucleon, which can 
be exactly compared to a^p computed here if DM equally couples to p and n. This is true for the example 
presented here where DM has universal scalar coupling to all quarks. There can be 0(1) difference if DM 
couples to neutrons and protons differently, but this does not alter the main point here. 
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B. Basic Mechanisms Affectine \hherm 

For some models, a dark matter candidate does not naturally satisfy {(t\v\) therm = 1 pb 
and requires some enhancement in order to yield a thermal relic abundance. This simulta- 
neously modifies the natural ratio in eq. (|22l) . The two basic mechanisms are: 



S-Channel Resonance : If dark matter annihilation is enhanced by an s-channel res- 
onance, eq. (!22l) can be satisfied because of kinematic suppression in the denominator of 
eq. fl26|) . This solution requires one to tune the mass relation between DM and mediator: 
2m OM ~ M. With highly tuned 2m£)M = M, resonance enhancement depends on the 
width-to-mass ratio, ~ M^. As an example consider that for a Z-resonance, this factor is 
about 10'^ and for a Higgs resonance {rrihiggs = 100 — 200GeV), this enhancement ranges 
from 10^ — lO^''. In practice, some models require tuning at the few percent level to satisfy 
eq. (122]). 

Coannihilation : If dark matter coannihilates with other states close to its mass, eq. ( !22|) 
can be satisfied because {cr\v\) therm will be dominated by coannihilation with its degenerate 
partner, rather than self annihilation. This also requires tuning mass relations in the model 
so that DM has partners with almost degenerate masses. It is most useful when DM self- 
annihilation is f ^ suppressed, while coannihilations are unsuppressed. Then this can increase 
{o'\v\)ann by about l/t>^ ~ 30. This is typically not enough to solve the O(IOO) tension we 
mentioned above, and therefore suppressions of asi are generally required. 

If a dark matter candidate does naturally have {cr\v\)therm = 1 pb, then one or a com- 
bination of the following mechanisms usually explains how eq. ( 122|) is satisfied. The basic 
principle underlying these mechanisms is the absence or suppression of the SI coupling to 
light quarks. Since eq. (!22l) arises from SI DiDt experiments, it is possible that the only 
operators governing annihilation, when crossed, give rise to SD scattering or do not couple 
to light quarks at all. While the crossed annihilation channel may not give rise to SI scat- 
tering, other operators typically do generate SI scattering. In these cases one still expects 
that {cr\v\) therm = 1 pb implies eq. ( I30l) if SI scattering is generated by other operators. Of 
course the leading operators that do mediate SI scattering should still obey DiDt bounds. 
Therefore, even if the crossed annihilation diagram does not mediate SI scattering, one of 
the following mechanisms is typically required to suppress asi- 
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Kinematic Suppression of SI Coupling: Kinematic differences for annihilation vs. 



scattering can suppress a si relative to {cr\v\) therm- Suppressions can arise for certain oper- 
ators from either the slow current-day velocity of dark matter {e^ = {^^^r^Y ~ 10^^) or the 
low momentum of quarks in the nucleon {eqcD = (j^^^ ~ 10"^). This can be in addition 
to the automatic (^)^ factor in eq. (^71) . which is always present. We give a brief summary 
and explanation of the kinematic suppressions in Appendix-A (Table-lV). More interestingly 
for certain operators, SD interactions have a large <7sd, even when asi is suppressed. Models 
with significant SD scattering might be probed with experiments in the near future, with 



exclusion limits as low as (Tsd ~ 10 ^pb, for instance, in ICE CUBE 



44|. 



Suppression from Flavor Dependent Couplings: If the DM coupling to Standard 
Model light quarks is suppressed so that asi is suppressed, then eq. ( 122|) is readily satisfied. 
Meanwhile, {(y\v\)therm remains unsuppressed in the presence of other efficient channels: 
heavy quarks, leptons, or W/Z. Such possibilities include: DM coupling via a Higgs-hke 
mediator which leads to Yukawa suppression or Z-boson mediation with mixing suppression. 
Yukawa suppression is a natural possibility that is incorporated in many known models. 
However, without incorporating additional suppression mechanisms, Yukawa suppression 
alone has 'limited' power, and predicts cxsi ~ 10^^ — 10^^ pb (ey ~ 10"'*) which will soon 
be probed by XENONIOO and XENONIT experiments. We can understand this as follows: 

Going back to eq. fl26l) and replacing g2 by suppresses the annihilation cross-section for 
the u-quarks. Taking Yukawa couplings (for all quarks) into eq.(l27ll (replace the universal 
92 by Vg): 

2 



47r (m^ + mJ2 \ ^ m/''''^^ ' ^ m/''27' 

q=u,d,s ^ q=c,b,t 



1 ml (^i)2 mp 2 



A. 

with a WIMP and Higgs mass ~ lOOX GeV and g2 = gEW, with eq. (!26|) . we get <7DiDt ~ 
^(^)^^(o"|f Dtfeerm ~ 10~''^cm^. So we see that this type of model is right around the 
reach of XENONIOO or XENONIT. 

Inelastic splitting: In this type of model, DM has a heavier 'excited' partner, which 
opens up the possibility of inelastic scattering. If dark matter's elastic scattering to itself is 
suppressed or vanishes eq. (122!) can be satisfied. In these cases, scattering mainly proceeds 
via inelastic scattering to a partner with heavier mass. Enough mass splitting between DM 
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and its heavier partner gives a kinematic barrier, i.e. the scattering can occur only when DM 
carries enough kinetic energy. Therefore asi is suppressed by the number density distribution 
of DM at high velocity, which can be exponentially suppressed if SI scattering requires a very 
high velocity. An example of this is adding a Majorana mass which splits Dirac particle/anti- 
particle states into two mass eigenstates such that SI scattering is dominated by inelastic 
scattering. For Majorana masses larger than 1 MeV this is sufficient to evade all DiDt 
bounds. With a mass splitting ~ lOOkeV, such models have been recently well explored for 
reconciling DAMA signal and null results from other DiDt experiments. 

Annihilation to Dark Sector States : It is possible that DM annihilates dominantly 
to non-Standard Model states and only couples to the Standard model via small mixing. In 
this case, the correlation between {<7\v\) therm can disappear completely. Such type of models 
)een recently well explored in light of explaining anomalies from PAMELA, FERMI 



have 
etc. 



45| . In these models DM annihilates dominantly to light dark sector fields which 



gives rise to right thermal relic density. Elastic scattering off nucleons is typically highly 
suppressed by small mixing between Higgs and light scalar unless proper mass splitting could 
accommodate a sizable inelastic scattering rate 45|, |46 | 

Non-Thermal Production : The correlation between DiDt and annihilation cross- 
section disappears when dark matter is not thermal relic. We will not consider further 
details of these models here, but note that this simply eliminates the bound eq. fl22|) . 
Despite the popularity of thermal WIMP models, there are many well-motivated non-WIMP 
candidates such as the axion or gravitino LSP. Many of these models are 'super-weakly' 
interacting both for DiDt and indirect detection. 



V. SOME SPECIFIC WIMP MODELS 

In Section IIIII we have restricted to our attention to various dark matter annihilation 
modes as if they occur 100% of the time into a particular final state with {cr\v\)ann = 1 pb. 
However, dark matter candidates typically can annihilate into various final states in different 
ratios in a concrete model. In this section we select a small subset of specific dark matter 
candidates and consider the extent to which they can be probed by the various antideuteron 
searches. Each of these models avoids the naive correlation between DiDt and thermal freeze- 
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out due to some combination of the mechanisms mentioned in the previous section. For each 
model, we hst the leading operators that govern the particle's present day annihilation, its 
SI scattering, and SD scattering. We briefly describe the basic features of each model and 
which of the previous mechanisms is responsible for that model being kept within the bounds 
of DiDt experiments. We'll see that several models that are difficult to detect directly can 
be detected through annihilation to antideuterons. 

The various cross-sections derived from dark matter/Standard Model operators may come 
with certain suppressions. For {a\v\)ann, "P-wave" (v^) suppression can easily be inferred 
from the C and P properties of the dark matter bilinear part of the operator or CP, J 
conservation of initial-final states. Since the expectation value of various quark bilinears 
in the nucleon state determine the behavior of the SI and SD scattering cross-sections, 
kinematic suppressions for these can be determined simply by looking at the fermion bilinear 
to which dark matter couples. In what follows in this section, we will only state whether or 
not the leading operator has kinematic suppression and mention what type of suppression 
it is (Yukawa: ey, velocity: Cy, or nucleon energy scale: egcn)- We provide a detailed 
discussion of these suppression factors and their origin in Appendix A. Also, Appendix A 
contains a description of the compact notation used here for listing dark matter/Standard 
Model interactions. 



A. The Models 

1. Focus-Point SUSY: Low Mass 



F.P (1) 


Xl^l^'xht^l^'t (ANN) 


XXff (SI) 


xri'xfri'f (SD) 



This SUSY dark matter candidate is mixed Higgsino and Bino having mostly Bino 
component. A large annihilation cross-section is achieved to top quark final states via 
Z-exchange (though W/Z final states are also relevant). This "S-wave" annihilation 
channel is helicity suppressed for light quarks because the initial and final states must 
be CP odd and have J — 0. However, in "S-wave," the Standard Model fermion 
interaction produces only the J = 1 (5" = 1) final state. A helicity flip (mass insertion) 
is required to conserve CP and total angular momentum. Note that for fighter DM 
masses the top mass does not greatly suppress the ampfitude. So this cross-section 
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determines the thermal rehc abundance. SI scattering is dominated by Higgs exchange 
and is Yukawa suppressed. SD scattering arises from the crossed annihilation channel 
(with light quarks) and is not suppressed. 

2. Focus-Point SUSY: High Mass 



F.P (2) 


X-f'''y^x{^WdW)+^ (ANN) 


xxff (SI) 


xri'xfri'f (SD) 



This SUSY dark matter candidate is also mixed Higgsino and Bino but with more 
Higgsino component. Because the Higgsinos couple W-bosons, a large annihilation 
cross-section is achieved to W/Z final states (the top quark final states are also rel- 
evant). The "S-wave" annihilation channel is unsuppressed and determines thermal 
relic abundance. SI scattering is dominated by Higgs exchange and is Yukawa sup- 
pressed. SD scattering arises from Z-boson exchange and is mildly suppressed due 
to a slight degeneracy between the higgsino mixture of the dark matter state j^^]. 
Essentially the SD coupling of the neutralino to the Z -boson is ~ — l^-ffdP 

where Zh is the fraction of LSP having Higgsino of a particular variety and the minus 
sign is from the generator of SU{2)l. As DM gets heavier, the Higgsino state that 
is the LSP is equal mixture of the two weyl Higgsinos and this couplings becomes 
suppressed. 

3. SUSY BINO: with Coannihilation 



Coann 


Xl^l^xbl^l^h (ANN) 


xxff (SI) 


xri'xfri'f (SD) 



This SUSY dark matter candidate is Bino like. The annihilation cross-section is to bb 
final states, is helicity or "P-wave" suppressed, and too small to yield the correct relic 
abundance. A thermal relic abundance is achieved via coannihilations with a mass 
degenerate species (usually stau) in the early universe. Higgs exchange dominates the 
SI scattering and is ey suppressed. SD scattering arises from Z-boson exchange. 

4. SUSY BINO: with s-channel Resonance ("A-funnel") 



A-funnel 


Xl'xh'f (ANN) 


xxff (SI) 


xr^'xfYl'f (SD) 
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This SUSY dark matter candidate is Bino like. A large "S-wave" annihilation cross- 
section is achieved from resonant s-channed CP-odd Higgs boson exchange assuming 
that m^jjvf ~ 2m^. This can yield the correct relic abundance. The crossed annihi- 
lation diagram mediates subdominant e^cy suppressed SD scattering. Dominant SD 
scattering for DiDt is through Z-exchange. SI scattering is via Higgs exchange, and is 
ey suppressed. 

5. UED Model: KK fi^^) 



UED 


{eVdVrj-fJ (ANN) 


{eVdVrj^J (SI) 


{eVdVnh.Y^f (SD) 



This UED [48] dark matter annihilates to fermions via t-channel KK fermion exchange. 
"S-wave" annihilation to quarks is unsuppressed and can yield the correct thermal relic 
abundance. Crossed versions of the annihilation diagrams yield SI and SD scattering 



49| . SI scattering via this operator^ is kinematically suppressed because only the 
time component of the dark matter momentum is unsuppressed, implying that only 
the f'Yf component of the fermion part of the operator is relevant and yields 
suppressed SI scattering. The SD scattering is not kinematically suppressed, however, 
since the mediator is a KK (1 TeV) state rather that a Z-boson. The SD scattering is 
roughly 10^^ of the typical SUSY SD scattering. 

6. Little Higgs with T-Parity Model 



LHTP 


{VV){WW) (ANN) 


VVff (SI) 


(eVdVrjj.j'f (SD) 



Little Higgs models with T-Parity (LHTP) 50l-l52| contains a spontaneously broken 
"mirror" gauge symmetry and the dark matter is a massive "mirror" photon of its 
f/(l) factor. In these models dark matter can "S-wave" annihilate through a Higgs 
resonance and produce the correct thermal relic abundance, but the natural anni- 
hilation cross section is typically too small and requires itldm = 2mHiggs in order 



^ We have listed only one of the operators one gets from decomposing the t-channel fermion exchange 
diagram into operators with dark matter and Standard Model billinear operators. However, the other 
operators have equal or more kinematic suppression. For example {VV)'^ {f^^dyf)- which is associated 
with eqcD, while cqcd ^ for a weak scale DM 
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enhance the s-channel annihilation to the correct value. This induces an s-channel 
resonance suppression of the ratio: isidlti^iiri _ Since the annihilation is not to fermions 



the SI and SD scattering are determined by different operators than annihilation 53 |. 
The SI interaction is dominated by Higgs exchange and therefore gets a Yukawa sup- 
pression in addition to the suppression relative to the s-channel resonance . The SD 
interaction is similar as in UED in addition to another factor of 10~^ suppression. 
This suppression arises from the fact that in Little Higgs Models with T-parity, the 
Standard Model is charged under a spontaneously broken gauge group with "mirror 
hypercharge" (Y' = since the dark matter is the massive photon of that spon- 
taneously broken symmetry, then it's SD elastic scattering cross section off quarks is 
suppressed by Y"^. 

7. Warped Extra-dimensional Model: KK LZP u'^^ 



LZP 


xrxilf.t (ANN) 


xrxhj (SI) 


xi^i'xh.i'f (SD) 



In models of warped extra dimensions it is possible to have dark matter that couples 



mainly to tops 5J, |55(| or more generally, a Dirac fermion with flavour dependent 



couplings to quarks via Z' gauge bosons |56|. In these cases, annihilation to It is 
unsuppressed and gives the correct relic abundance. SI scattering is not kinematically 
suppressed, but is small due to a small Z/Z' mixing. SD scattering is suppressed in 
the same way. This suppression is similar in principal to "Yukawa" suppression where 
flavor dependent couplings play an important role. 

8. Neutral Singlet Scalar Mixing with Higgs 



Scalar 


0t0/i/i (ANN) 


0V/7 (SI) 


(SD) 



One of the simplest dark matter mode 



interactions with the Higgs boson 57 



. s of all is a neutral stable scalar having quartic 
SoJ. For mass above m^, DM annihilates to 
Higgs bosons, resulting in the correct thermal relic abundance. W-bosons can be an 
important final state for some scenarios. The SI scattering is Yukawa suppressed and 
there is no SD scattering. 

9. Doublet/Singlet Scalar Inelastic Model 
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Doublet /Scalar 


{<P^d^r_h,f (ANN) 


(SI) 


{^^d^r_h,j'f (SD) 



A light scalar doublet /singlet mixture whose mass is just below the W-mass threshold 
can provide a correct relic abundance. This model is a variant of models appearing in 
60| . We will leave a discussion of the details of this model for Appendix B and only 
mention the essential features here. The dominant annihilation in the early universe is 
to quarks and this is "P-wave" suppressed, but for light enough dark matter can yield 
the correct relic abundance. Present day annihilation is dominated by annihilation 
to two (slightly) off-shell W-bosons so that {(T\v\)ann = 0.1 pb. SI elastic scattering 
is absent, but inelastic scattering is possible. This feature can be used to explain 
the apparent conflict between DAMA annual modulation signal and null discovery in 
other DiDt experiments 611]. General strong neutrino bounds on inelastic DM from 
solar capture ( 66|) are somewhat relieved in this model as a result of lighter mass 
and mixing with singlet. In this particular model the SD scattering is mediated by 
Z-boson exchange and is of typical size. 



10. SUSY Wino: non-thermal 



Wino 


Xl''l^x{^WdW)+^ (ANN) 


xxff (SI) 


xri'xfri'f (SD) 



lave present day cross sec- 



67l |. It has a large present 



Finally, we consider one dark matter candidate that does not 
tion equal to that of a thermal relic, a Wino-like SUSY LSP 
day annihilation cross-section. SI scattering arises from Higgs exchange and is Yu 
pressed [68|. SD coupling arises from Z-boson exchange and is also suppressed 



cawa sup- 



69|. This 

model predicts a very large antideuteron signal for all experiments mainly due to its large 
annihilation cross-sections (10^ pb). It is plausible that propagation parameters other than 
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7l| and can 



those quoted in this paper can be used for the astrophysical background 
be in agreement with anti-proton backgrounds. While relaxing the allowed propagation of 
background antiprotons one does not significantly reduce the propagation of dark matter 
antiprotons or antideuterons. Thus a Wino candidate, such as this one, should be robustly 
detected or ruled out by future antideuteron searches. 
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B. Model Summary 



The leading present day annihilation channels for all of these models are not "P-wave" 
suppressed and have significant annihilation to final states that produce hadrons. Therefore 
they will contribute to the population of antideuteron cosmic rays. Picking particular pa- 
rameters for each model, we can use our previous general analysis to determine each model's 
dark matter contribution to the antideuteron cosmic ray spectrum and calculate an expected 
number of antideuterons for the GAPS and AMS-02 experiment. Our results are summa- 
rized in Table lllli We list the mass and present day annihilation cross-section along with 



the fraction of annihilations to final states yielding hadrons. 



or the supersymmetric models 



we have used the numerical package DARKSUSY 5.0.4 |62l . |63| to check that the choice of 
parameters yields the correct relic abundance. The choices of parameters are similar to the 
parameter space described in 6J] for the Focus Point regions of parameter space and 65 1 
for the Coannihilation and "A-funnel" regions of parameter space. We list the predicted 
number of antideuteron events for the GAPS (ULDB) experiment as well as an estimate for 
a satellite mission, GAPS (SAT). The number of events predicted should be compared to 
the Ncrit column of Table [TTl in order to determine whether or not the predicted signal is 
significant. Finally, we list the expected order of magnitude for SI/SD elastic scattering off 
of nucleons, which reflects the kinematic and coupling suppressions mentioned above. 

We have not included the AMS-02 and GAPS(LDB) experiments in this list because 
they are not sensitive enough to detect most of the models for the parameters chosen here. 
However, Wino dark matter is a notable exception. For m^jM = 200 GeV, this predicts 10 
(12) events in the high (low) energy AMS-02 experiment and 8 events in the GAPS (LDB) 
experiment. Similarly, for moM = 300 GeV, this predicts 2 (1) events in the high (low) 
energy AMS-02 experiment, and 1 event in the GAPS (LDB) experiment. While AMS-02 
and GAPS (LDB) are not promising for probing the particular thermal relics discussed in 
Table IIIIl based on our model independent analysis in Section IIIIl we can see that the mass 
reach for AMS-02 is 150 GeV(110 GeV) for the low (high) energy window. Also, while 
GAPS(ULDB) will have greater sensitivity than AMS-02, AMS-02 will remain the only 
probe of antideuterons with kinetic energies Tg, > 1 GeV. 
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Model 


TTlDM 

(GeV) 


0" f l/pb 










Nsa = 2 

N5a = 5 

(ULDB) 


N3a = 10 
N^a = 16 

(SAT) 


asi/n/ph 


asD/n/ph 


SUSY F.P (1) 


190 


0.67 


0.2 


0.02 


0.73 





4 


50 


10-8 


10-4 


SUSY F.P (2) 


772 


0.33 


0.55 





0.38 








4 


10-8 


10-5 


SUSY coann 


148 


0.17 





1 








1 


14 


10-8 


10-6 


SUSY A-funnel 


163 


0.6 





0.92 








2 


33 


10-8 


10-6 


UED 5(1) 


900 


0.6 





0.19 


0.16 


0.02 





3 


10-8 


10-6 


UED B(i) coann. 


600 


0.6 





0.19 


0.16 


0.02 





4 


10-8 


10-6 


LHTP 


200 


0.8 


1 














12 


10-12 _ 10-10 


10-10 


LZP 


300 


1 


0.06 





0.94 





3 


41 


10-9 _ 10-6 


10"^ - 10-4 


Singlet (scalar) 


200 


1 











1 


2 


36 


10-8 





Doublet /Singlet 


75 


0.1 


1 











3 


49 





10-4 


Wino (non-therm) 


200 


65 


1 











50 


730 


10-12 - 10-8 


< 10-11 


Wino (non-therm) 


300 


30 


1 











5 


75 


10-12 _ 10-8 


< 10-11 



TABLE III: Predicted number of antideuterons detected in various experiments for different dark 
matter models. Here where "/" is some final state, signifies the fraction dark matter annihila- 
tions that yield this particular final state. Also, stands for all light quarks (including b-quarks) 
and includes only top-quarks. Otherwise these are the same parameters as in Table HIl We also 
list the typical range for SI and SD elastic scattering off nuclei. 



VI. CONCLUSIONS 

In this paper, we have considered the mass reach of the antideuteron cosmic ray search 
experiments AMS-02 and GAPS. The basic point can be summarized by stating that for 
typical thermal dark matter annihilating with a present day annihilation cross-section of 
{(^\v\)ann = 1 pb, GAPS (ULDB) expects to detect antideuterons of dark matter origin at 
greater than 99.7% confidence for masses up to m^Af = 250 GeV in the case of annihilation 
to h^h^, niDM = 400 GeV for tt, tudm = 250 GeV for qq, tudm = 160 GeV for WW, 
and rriDM = 500 GeV for gg final states, assuming MED propagation parameters. At 
5cr(~ 99.99%C.L.) GAPS (ULDB) expects to discover antideuterons of dark matter origin for 
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DM masses up to mciM = 150 GeV for mcM = 220 GeV for tt, m£,M = 150 GeV for gg, 
^DM = 120 GeV for WW, and m^iM = 300 GeV for gg final states. Therefore antideuterons 
are a significantly more sensitive way to probe the hadronic annihilation channels of dark 
matter than antiprotons. A GAPS experiment on a satellite, obtaining an order of magnitude 
improvement in sensitivity, would provide the same mass reach at 5a significance that GAPS 
(ULDB) achieves for 2a significance and could see primary antideuterons for dark matter 
masses of 1 TeV. GAPS (ULDB) is virtually background free, and satellite level sensitivities 
are expected to see only a few background events. AMS-02 and the GAPS (LDB) are 
also sensitive to antideuteron cosmic rays of dark matter, and the mass reach is roughly 
moM = lOOGeV. 

We have computed injection spectra according to the non-factorized coalescence model 



261] which more accurately takes into account angular correlations in antideuteron forma- 



tion when hadronization is modeled with PYTHIA. Since this has, correctly, increased the 
prediction in antideuterons at higher energies we have produced flux results that are some- 
what flatter at higher energies than previous analysis. While this does not make a great 
difference for low energy antideuteron searches, it does affect the higher energy AMS-02 
detection rates. For a fixed annihilation rate of {a\v\)ann = 1 pb and MED parameters, we 
find a 100 GeV dark matter candidate predicts 1 event for annihilation to bb. Since dark 
matter can be as light as 50 GeV without violating the antiproton bound, this high energy 
band may still be relevant for dark matter searches. 

The dominant uncertainty comes from uncertainties associated with cosmic ray propaga- 
tion from dark matter annihilation . The minimum flux can be an order of magnitude less 
than the flux inferred from using the MED parameters used to make the plots in this paper. 
This results in degrading the mass reach by a factor of 3. However, the mass reach, while 
less, remains an order of magnitude better than the bound from antiprotons. Because the 
model of antideuteron production is fit to the ALEPH measurement of Z-boson decay to an- 
tideuterons, we believe the hadronization uncertainties are 0{1), which are small compared 
to propagation errors. 

We have also described the general tension between dark matter as a thermal relic and 
bounds set by Direct Detection experiments. We have listed the basic mechanisms that are 
at work in many different models of dark matter responsible for eliminating this superficial 
conflict. These mechanisms typically manifest themselves as different suppression factors for 
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SI scattering. In Appendix A we tabulated these suppression factors associated with different 
Standard Model fermion bilinears, when their expectation value is taken in a nucleon state. 
We then used Appendix A to understand different dark matter models by their leading dark 
matter/Standard Model operators which determine annihilation and elastic scattering off 
nucleons. We identified which mechanisms/suppression factors are present in each model 
alleviating the conflict between the dark matter as a thermal relic and Direct Detection 
bounds. Finally, for each of these models, we have picked a particular dark matter mass 
(and model parameters) that satisfies Qh'^ — 0.12 and determined the antideuteron signal 
predicted in the GAPS (ULDB) and GAPS (SAT) experiment. Our selection of models is 
meant to be representative of typical models but by no means complete. We expect that 
GAPS (ULDB) will provide a first real test for hadronic annihilation channels for dark 
matter with masses smaller than 100 GeV. For GAPS (SAT), the possibility exists to probe 
hadronic annihilation for masses up to itidm — 1 TeV. 
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Appendix A: Operator Properties Relevant for Dark Matter Detection 

As commented in section IV, a general DM candidate's annihilation cross-section, SI 
scattering, and SD scattering depend on the operators involved in interactions with the 
Standard Model (e.g. whether there is suppression in SI DiDt rate or IdDt rate). In this 
Appendix we shall classify dark matter scalar, fermionic (Dirac and Majorana), and vector 
boson bilinears according to their transformation properties under Charge Conjugation (C) 
and Parity (P) and derive kinematic suppressions for the expectation value of bilinear opera- 
tors between the vacuum and a non-relativistic (^ <S 1) particle anti-particle state. This will 
be relevant for determining whether "S-wave," "P-wave," or helicity suppression plays a role 
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in dark matter annihilation. We present this for completeness; much of what appears here 



can be found in standard texts 
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2|-|74||. We will also tabulate kinematic suppression factors 
for quark bilinears in a nucleon state, which will be important for interpreting bounds from 
Direct Detection (DiDt) experiments. 

We list a basis of bilinear operators for fermions in Table \IV\ for scalars in Table 
IV] and for gauge bosons in Table |VT1 These tables are to be used in Section |V] to 
infer the kinematic suppressions for annihilation or SI/SD scattering for particular dark 
matter/Standard Model operators. These tabulated bilinear operators can be interpreted 
as being dark matter or Standard Model field bilinears. All possible 4 point interactions 
between Standard Model and dark matter can be written as a product of two of the 
bilinear operators tabulated here in form of OdmOsm- Any dark matter/Standard Model 
operator can be brought into this form. For example, {xio- + i>lb) f){f {o- — bl5)x) for SUSY 
neutralino-SM fermion interaction via t-channel exchanging sfermion, can be reorganized 
using Fierz identities to a combination of operators of the form described above. Restricting 
to CP-conserving interactions simplifies the number of 4 point interactions, however the 
number of operators is still quite large. We restrict our attention to interaction operators 
up to dimension-6. Some short-hand notations are used in order to fit into a compact 
table, in particular for vector boson fields, which will be explained. One common notation 
in form of (x^x)± is short-hand for two bilinears with definite C parity: x^X + X^X^ 
i{xAx — X^Xji respectively. Obviously for a real field only the + one survives. Another 
= where (-1)° = +1, (-1)^ = = 1, 2, 3). 



In the case that bilinears are of fields describing dark matter, it is a simple exercise to 
establish whether or not present day annihilation is "P-wave" suppressed or not based on C 
and P properties of the bilinear operator and and "S-wave" state (at the end of this appendix 
we comment on some particular cases where CP, J conservation between initial and final 
state needs to be considered to determine "P-wave" suppression or helicity suppression at 
"S-wave" ) . 

Consider the massive particle/anti-particle state: ^i, Si; ^i, -52). This state can be 
decomposed into a superposition of orbital angular momentum eigenstates. Since dark mat- 
ter annihilates in the non-relativistic limit, we will neglect all but the / = 0, 1 states, also 
called "S-wave" and "P-wave," when the states are initial states of a 2 — > 2 annihilation 
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process. For "S-wave" or "P-wave" states, the state then has well defined C and P trans- 
formation properties given by: 



C = (-1)^+^,P = (-1)^ (Scalars/GaugeBosons) 
C = (-1)^+'^,P = (-1)^+^ (Fermions) 



(Al) 



(A2) 



(A3) 



Where "L" is the orbital angular momentum quantum number and "S" is the total spin 
quantum number. Since dark matter annihilation to the standard model will always involve 
expectation values of the form: 



where O is some dark matter bilinear operator. Then we may immediately use the C and 
P properties of the state and the operator to determine whether these expectation values 
are zero or not for an "S-wave" initial state; in order to be non-zero the quantum numbers 
of the state and operator should match. 

For example, the fermion bilinear xx always contributes "P-wave" suppressed annihilation 
because there does not exist a fermion/anti-fermion state that is C and CP even with zero 
orbital angular momentum. In the tables " /" indicates that eq. flA4p is non-zero for an 
"S-wave" state and a given dark matter billinear. Furthermore, it is not uncommon for dark 
matter candidates to be their own anti-particles. In this case all two particle dark matter 
states are even under Charge Conjugation and only the operators with "-I-" subscripts yield 
no n- vanishing contributions to either annihilation or SI/SD scattering. 

If one considers the case that the fermions in the fermion billinears of Table IIVI are 
Standard Model light quarks, then the expectation value of these operators in a nucleus 
state determines important features for ctsi/sd'- whether it is SD or SI and whether it has 
certain type of suppression as indicated by ey, ev, cqcd as explained below (ey, also apply 
when the fermion in the bilinear is DM that scatters off nucleus.). In the first two rows ("SI', 
"SD" ) of Table HV] we summarize the contribution of the operator for momentum transfer 
(g) which for DiDt is typically g oc f ~ O(keV) ^ mN,mDM- / means no suppression; 
indicates strictly vanishing; indicates a suppression ~ ti^ ~ 10~^ for aDiou ^qcd indicates 
a suppression ~ (■^7^)^(~ 10"^ for rriDM ~ lOOGeV); ey indicates that a possible Yukawa 



{pi,h,si;p2,h,S2\O{x,x)\0) 



(A4) 
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SI 
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ev 




ev 


SD 







ev 
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ev 


^QCD 
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+ 


+ 




+ 






=F 


± 


p 


+ 




(-r 


-i-r 






i-r" 




s-wave 





/ 


/ 






/ 


+ :/,-: 


+ : 0, - : / 



TABLE IV: Here is a list of fermion billinears, up to dimension four, having definite C and P 
transformation properties. Listed is each operator's suppression factor when its expectation value 
is taken in a nucleon state. Also listed is whether or not the operator allows "S-wave" annihilation 
for a dark matter state. 





0t0 






c 


+ 


± 


± 


p 


+ 






s-wave 


/ 


+ :/,-: 


+ :/,-: 



TABLE V: Here is a list of scalar billinears, up to dimension four, having definite C and P trans- 
formation properties. Listed is whether or not the operator allows "S-wave" annihilation for a dark 
matter state. 



suppression ~ (^j^)^ ~ 10~^ if the operator originates from integrating out a Higgs-like 
mediator. It can be seen that typically either SI or SD has suppression, but not both, 
which results from the fact that in the q — > limit either or cr* part of matrix element is 



picked up. Details can be found in [39|, l75 |. 



For vector boson billinears we adopt the following notation: 

VV = V+V-, iyVY'' = V^^V-\ {eVVy' = e^'^'P^V^V- , (A5) 
{VdVy = V+^d^V.^, V+'^d^'V.^, V+'^d^V-", {eVdVY = e>"'f"'V+^dpV^^, 
{VddV)^"' = V+^'d'OPV-, V^^'d^V-", V^^dpOPV-", V^f'dpd^'V-'' , 

Vd'^V = IZ+'^SV", eVddV = e^^P^V^d^dpY- 
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[eVVf^ 


{VdVf^ 






{eVddVf^ 




{eVddV)± 
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p 


+ 
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+ 




s-wave 


/ 


/ 


/ 


+ :/,-: 



TABLE VI: Here is a list of Massive Spin-1 billinears, up to dimension four, having definite C and 
P transformation properties. Listed is whether or not the operator allows "S-wave" annihilation 
for a dark matter state. 



A dark matter /Standard Model operator of the form OdmOsm is in general a combina- 
tion of two of the bilinears listed above. Even if one restricts to CP conserving operators, 
the list of possible operators up to dimension six is well above 50. Therefore we do not list 
these operators here. However, given such an operator, kinematic suppressions for DiDt can 
be read off directly by looking at which SM fermion bilinear appears in the operator and 
reading off its suppression factor from Table IIV[ Similarly, in order to determine whether 
or not annihilations are "P-wave" suppressed, one typically only needs to determine if the 
cooresponding dark matter bilinear has has non-zero expectation value in a two particle 
"S-wave" dark matter state. However, in some cases this information in combination with 
selection rules from CP and J conservation is required in order determine whether or not he- 
licity suppression is present. We give two brief examples here. First, in order to understand 
"S-wave" annihilation via the operator Xl^l^ xf 1^ {l^) f {^"^ important operator for SUSY 
Bino LSP annihilation), notice that xi^l^X has C = +1 and fermion pair state has P = — 1 
and determines the initial state should beatCP = — 1,S' = state at L = 0. To conserve 
J = final state needs to have opposite spin orientation. But for massless fermion ampli- 
tude for this helicity structure vanishes for state connected to /7^(7^)/. Therefore a mass 
insertion is needed to flip helicity, which implies a Yukawa-like suppression for annihilation 
to light SM fermions. Another example is, Xlnl^xi^dVY- Again for the same reason as in 
the first example, at L = initial fermion state should be at J = 0, C = +1, P = —1, which 
selects fi = component of X7^7^X to satisfy P = —1. Accordingly /i = or "S-wave" 
component is selected for {VdV)':!L, which vanishes as summarized in the table. Therefore 
the initial state must be in a "P-wave" state and (cr|f is suppressed. 
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Appendix B: Doublet /Singlet Scalar as Inelastic Dark Matter 



Motivation 

Inelastic dark matter (IDM) remains one of the most promising models to reconcile 
DAMA and other DM direct detection experiments. But since a relatively large nucleon 
scattering (a > 2 x 10^'^'^cm^) is needed to explain DAMA signal, it recently received 
stringent bound considering neutrinos from DM capture/annhilation in the sun [g^. The 
cross-section per nucleon times branching to WW is constrained to be a < (9(10~^^cm^). 
The constraint is particularly severe with the pure SU{2)l doublet IDM model considered in 



601 1 ■ Combining this with the thermal relic density requirement, we need a ~ 7x 10^'^^cm^ to 
fit all data. This shows O(IO^) tension with neutrino bounds. For this reason it is important 
to ask if IDM may be realized in a different framework. 

In this Appendix we shall propose a different IDM candidate than the scalar mentioned 
above. We consider a light scalar Dark matter with M^m < that is a doublet/singlet 
mixture. Our aim is to show that this can be compatible with all the above constraints 
simultaneously. Futhermore, as described in the the main text, this IDM scenario favors low 
energy antideuteron production, which is why we bring such detailed attention to it here. 
There are several motivations for considering DM lighter than my/. First, according to 



661], in the case where once two onshell W annihilation channel is closed, the constraint 
on a is relaxed to be 0(10"^*^ cm^) (unless r dominates annihilation). Meanwhile lighter 
mass implies larger DM flux which allows fitting DAMA data with lower a. So the tension 
between DAMA IDM and solar neutrino bound can be alleviated from two directions. For 
the case where DM is not much lighter than my/ the dominant channel is WW* producing 
3-body final states (as will be seen later, this is indeed the feature of our model presented 
here), though we will consider WW* are almost produced at rest. 

Light IDM as a pure SU(2) doublet underproduces dark matter from freezing out of 
thermal equilibrium. It was already found that doublet scalar (fermion) with mass of ~ 
500 GeV(l TeV) gives right thermal abundance. Roughly, the annihilation cross section 
scales This monotonic relation implies that it is impossible to get the right 



^ Not that when moM < mz/^, the scahng becomes Tn^u^Jm^. But for a pure doublet such low mass 
would hit collider constraint. 
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relic density for low masses. However an intriguing feature of a light scalar doublet is that 
the dominant annihilation channel at freezeout is to light SM fermion pairs via s-channel 
Z-exchange which has additional "P-wave" (f ^) suppression. This could compensate the 
increase effect on the early time {(y\v\) therm arising from lighter masses. From the operator 
analysis in Appendix A: the operator governing this is {<f)^d'^(j))^f'~f^{'y:^)f which "P-wave" 
suppressed. But as we will demonstrate by explicit calculation since ~ 1/20 at freezeout 
it still cannot provide a sufficient compensation for mw < 80GeV. A simple way to fix this 
is to mix the doublet with a singlet scalar. The necessary mixing angle can be determined 
by computing relic density of light pure doublet. 

Present day {cr\v\)ann is determined by the off-shell WW* final state because the phase 
space suppression is not as large as the velocity suppression for other channels. When the 
DM mass is lighter but close to niw, the phase space suppression is not severe so that the 
3-body annihilation via WW* can be dominant channel today with 0(0.01 — 1) suppression 
compared to {cr\v\) therm at freezeout. We will explicitly compute {<j\v\)ann of this 3-body 
annihilation and find that moM = 70 — 80 GeV works well for this purpose, and we shall 
take rriDM ~ 75GeV as our benchmark point, which gives (cr|f |) 

ann ~ 0. 1 ((T |f |)(^e^,^. It is 

worth mentioning that f ^ suppression for all fermion channels also efficiently suppresses the 
BR to r which could potentially impose stronger solar neutrino bound for IDM . 



The Model 



First it is instructive to compute the dominant {a\v\) therm at freezeout for a light pure 
doublet scalar which is annihilating to light SM fermions via Z-exchange 



{a\v\) 



1 



therm 



IGvr \ {s -ml + im^T^ 



Is - Ami r4(s + 2m'l){s - Ami] 



f 



7 

4ml 



3s 



(Bl) 



where for doublet Aw 



92 



2 cos 9w ' 



\f = Tf —Qf sin^ 6w The sum is over all light SM fermions. 



To a good approximation we set all mj 



0, r = 0, which leads to 



{a\v\) 



1 



therm 



167r |(s - m| + im^T^)\ 



We then estimate relic density using |76[|: 



Ami 



'A{s + 2m)){s - Ami] 



3s 



Y-^ = 0.26Agl/^mpm^3b/{x}) 
n^h^ ^ 2.82 X W^Y^{m^/GeV) 



(B2) 



(B3) 



37 



where b is the "P-wave" component from {cr\v\) therm, which is the coefficient of v"^ in eq. (lB2p 
after expanding for small f^. For the mass range we are interested in, it is reasonable to 
take h = 0.7, xj = 27, (7* = 9. We for DM of masses of 70 — 80 GeV, this pure doublet gives 
Qh^ = 0.01. Therefore as discussed earlier, to make critical thermal relic density we need 
to turn on mixing with singlet to reduce cross section by a factor of ~ 10. This implies 
that converting this model to a thermal relic requires mixing in a singlet candidate is mostly 
singlet with mixing angle sin^ ~ (0.05 — 0.1)^/"^ ~ 0.5. We will now construct such a model. 

Consider a SM singlet scalar S mixing with an SU{2)l doublet D with hypercharge 
Y = —1 via SM Higgs. We assume the inelastic splitting between real/irnaginary components 



of the scalar is explained by other physics mechanism (for instance in [60||) which decouples 
from our discussion here. We require a U{1) dm under which SM is uncharged, while S, D are 
charged to forbid large inelastic splitting from term like S'^h'^. This combines with proper 
mass inputs ensure stability of S despite of the vertex DhS*. For simplicity we set input 
coefficient for possible term |S'p/i^ to be at tree level. Radiatively generated terms of this 
kind are expected to have negligible effects by combining loop factors and Higgs Yukawa 
suppressions. The Lagrangian is: 

£ = \d^S\^ - ml\S\^ + \V^D\^ - m^|5|2 + {aDhS* + h.c.) (B4) 



where a is a dim-1 parameter (such mixing is similar to what we used in 60|). The mass 
matrix to be diagonalized: 



/ 7X1% av 


II 


1 


\ av mjj j 







S* D]\ ^ 1 (B5) 



Mass eigenstates can be written as {S will be our DM candidate, 6: mixing angle) 

S = cos^^-sin^D (B6) 
D = sin ^5 + cos (B7) 

Analytic solutions for mass eigenvalues and mixing angle are: 

1 



mf = -(m|) + m| — y m|) — 2m|,m| + m| + 4a^f ^) (B8) 
^2 — 2 ^"^^ ~'~ ~'~ \/ ~ 2m|,m| + m| + Aa^v"^) (B9) 



tan6' = ^^(m|, - m| + \J m^^ - 2m|)m| + m| + Aa^v'^) (BIO) 
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We require: 

1. S he DM candidate with mi < 80GeV 

2. We require sin 6 ~ 0.5 to get the right rehc density 

4. We require all input parameters to be Electro-weak scale 

A benchmark point that satisfies all the above requirements is: a = 15GeV, ms = 
85GeV, mD = lOOGeV, this input set gives 
rriDM = mi^ 75GeV, m2 ^ lOSGeV, sin6' ^ 0.5 



We have checked the compatibility of this dark matter candidate as IDM which explains 
DAMA data and other DiDt bounds along the lines of 60|, and find it fits at 99%CL with 
0" ~ 3 X 10~^°cm^. Such scattering rate is also more compatible with the solar neutrino 
bound which is ~ 10~'^°cm^ once highly energetic final state W's are forbidden and the r 
channel has negligible annihilation fraction. Meanwhile, as mentioned earlier, we find this 
benchmark point gives {cr\v\) therm ~ O.lpb today with WW* final states. 
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